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ABSTRACT 

We present parallax observations and a detailed model atmosphere analysis of 54 cool 
and ultracool (T e s < 4000 K) white dwarfs (WDs) in the solar neighbourhood. For 
the first time, a large number of cool and ultracool WDs have distance and tangential 
velocities measurements available. Our targets have distances ranging from 21 pc to 
> 100 pc, and include five stars within 30 pc. Contrary to expectations, all but two of 
them have tangential velocities smaller than 150 km s -1 thus suggesting Galactic disc 
membership. The oldest WDs in this sample have WD cooling ages of 10 Gyr, providing 
a firm lower limit to the age of the thick disc population. Many of our targets have 
uncharacteristically large radii, indicating that they are low-mass WDs. It appears 
that we have detected the brighter population of cool and ultracool WDs near the 
Sun. The fainter population of ultracool CO-core WDs remain to be discovered in 
large numbers. The Large Synoptic Survey Telescope should find these elusive, more 
massive ultracool WDs in the solar neighbourhood. 

Key words: techniques: photometric - stars: atmospheres - stars: evolution - white 
dwarfs - Galaxy: disc 


1 INTRODUCTION 


Given the finite age of the Universe, the first asymptotic gi¬ 
ant branch stars that for med now li ve as lo g (LfLpf) = —4.5 
white dwarfs (W Ds; iMesteJ 1_9oT_ Iben fc Tutukovl 
19841: IW ineet et al.l Il98 7l: iLiebert. Dahn fc Monetl 19881 : 


Fontaine. Brassard fc Bergeronl l20()fl l. Such WDs have 
temperatures below 4000 K (hence classified as ultracool) 
and they have been observed in deep Hubble Space Tele¬ 
scope ( HS T ) images of the halo glob ular clusters M4 and 
NGC 6397 (lHansen et al.ll200l . 120071 ). The oldest WDs in 
these two clusters are sal 1.5 Gyr old. 

Large-scale survey s such as the Sl oan Digital 


Sky Survey (SPSS; iGates et al. 120041: H arris et al.l 
20061. 120081: iKilic et al.l 120061. l201(il : IVidrih et al.l 120071 : 
Hall et al.l (2008k the UKIRT Infr a red Deep Skv Sur¬ 


vey (UKIDSS; Leggett et al.l l201ll : ICatalan et al. _ 20121 : 
iTremblav et al.l 12014 1 and SuperCOSMOS ( Hamblv et al.l 


* This work is based on observations obtained at the MDM Ob¬ 
servatory, operated by Dartmouth College, Columbia University, 
Ohio State University, Ohio University, and the University of 
Michigan. 

f E-mail: alexg@nhn.ou.edu 


Il999l : iRowell. Kilic fe Hamblvi l2008l l have identified the 
analogues of these ultracool WDs in the field. Since these 
field WDs are relatively bright compared to the globular 
cluster WDs, optical and infrared photometry in several 
bands can be easily obtained from ground-based telescopes, 
enabling us to model their spectral energy distributions 
(SEDs) accurately. This is important for understanding the 
different opacity sources in these stars, deriving reliable 
temperatures and ages, and also calibrating the faint WD 
sequences of globular clusters that usually rely on two filter 
photometry. 

The spectra of hydrogen-rich cool and ultracool WDs 
differ from those of their warmer counterparts because 
they show t he effects of the red-wing of the Lya opacity 
in the blue jKowalski fe Saumonl l2006l l and the collision- 
induced absorption (CI A) due t o molecular hydrogen in 
the near-infrared (lHanse n] I1999H . The latter shifts the 
peak of the SEDs of ultracool WDs back to the opti¬ 
cal wavelengths. Unfortunately, there are only three ultra¬ 
cool WDs in the field with parallax measurements. These 
are WD 0346+24(0, SDSS J110217.48+411315.4 (here- 


1 We note that this object is also known as WD 0343+247. 
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after J1102; Kiljc_et_aD 2012|, a nd re ferences therein) and 
LHS 3250 ( Bergeron fc Leggett) |200^ ~). The first two stars 
have SEDs that peak near 1 fim. On the other hand, 
the LHS 3250 SED peaks at 0.6 /im, representing an ex- 
treme case of CIA flux de ficit in the optical and infrared. 
iBergeron fc Leggett; ; 21II12 ' performed a detailed model at¬ 
mosphere analysis of LHS 3250 and demonstrated that 
LHS 3250 has a helium-rich composition, it is overlumi- 
nous, and undermassive. The best-fitting model and the 
parallax measurement ind icate a mass of only 0.23 M@ 
(IBergeron fe Leggettl 120021 '). This is somewhat problem¬ 
atic as all previously known low-mass WDs are DAs with 

hyd rogen-rich atmosph eres. _ _ 

iGates et alJ (120041 ') and lHarris et al.l (120081 ') as well as 
several other groups have identified about a dozen stars with 
SEDs similar to LHS 3250. In this paper, we present par¬ 
allax measurements and a model atmosphere analysis of 54 
cool WDs, including half a dozen ultracool WDs and several 
other cool WDs with significant infrared flux deficits. Our 
target s were selected from the cool and ultracool WD sam¬ 
ples o f lGates et ah d2004h. Vidrih et al.l d2007t) , lHarris et al.l 
(120081 1 and iKilic et al.l I 201ff ). and are biased towards WDs 
with significant infrared flux deficits. Parallax measurements 
allow us to accurately determine the distances, masses and 
consequently the cooling ages for these stars. Section 2 out¬ 
lines our observations including a description of our Bayesian 
approach to estimating distances. Section 3 describes the 
models used in our analyses followed by our results in Sec¬ 
tion 4. In Section 5, we discuss the ages and membership of 
the WDs in our sample as well as the implications of our 
results towards our understanding of WD evolution and we 
conclude in Section 6. 


2 OBSERVATIONS 
2.1 Parallax 

All our parallax data are from the 2.4m Hiltner telescope 
at Michigan-Dartmouth-MIT (MDM) Observatory on Kitt 
Peak, Arizona. We used a thinned SITe CCD (named 
‘echello’); at the /7.5 focus, each 24 /rm pixel subtended 
0.275 arcsec, giving a field of view 9.4 arcmin 2 . For all our 
parallax data, we used a 4-inch-square Kron-Cousins /-band 
filter, which did not vignette the CCD. Exposure times var¬ 
ied with the brightness of the object, but were typically a few 
hundred seconds. Our data were taken on numerous observ¬ 
ing runs between 2007 and 2011. Table [T] gives the epochs 
that each star was observed, and the number of exposures 
at each epoch. 

Our reduction and analysis procedures d i ffered only 
slightly from those described bv iThorstensenl (120031 ') and 
iThorstensen. Lepine fe Sharal d2008[) . As in the previous 
work, we corrected our raw parallaxes to absolute using 
colour-based distance estimates for the reference stars, and 
estimated uncertainties using the formal errors of the fit 
and the scatter of the references stars. In order to correct 
for differential colour refraction (DCR), we need to know 
the colour of both the programme star and the reference 
stars. In previous work we measured the colours, but for 
this work we used SDSS g — i colours and adjus ted th e DCR 
correction factor slightly to account for this. IThorstensenl 



Ka, sdssI ( mas y r_1 ) 



Figure 1 . Comparison between the proper motions measured 
at MDM Observato ry and those from the SDSS+USNO-B cat¬ 
alogue dMunn et ld] l2004l) for 42 of the 54 WDs in the current 
sample. We compare the absolute value of the proper motion in 
right ascension (|/traIi top) and in declination (|/'dj-;c|. bottom). 
The dotted line represents the 1:1 correlation. The dashed lines 
represent the ± 10 mas yr" 1 range. 


(200^) describes a Bayesian procedure used to estimate dis¬ 
tances from the available data, which combines the paral¬ 
lax measurement with an assumed space velocity distribu¬ 
tion and absolute magnitude range. We used a similar ap¬ 
proach here, but modified the prior information to be ap¬ 
propriate to the present sample. For the velocities, we used 
a composite distribution consisting of 60 per cent thin disc 
with ( U , V, W) = (39, 20, 20) km s -1 , 30 per cent thi ck disc 
with (U, V, W) = (46, 50, 35) km s -1 dChiba fc Beersll2000 h. 
and a 10 per cent probability of a still larger dispersion 
( U,V,W ) = (100,75,50) km s _1 . The absolute magnitudes 
of these WDs are likely to be in the range 11-18, so the abso¬ 
lute magnitudes were assumed to be drawn from a Gaussian 
centred on M B = 15 with a standard deviation of 4 mag. 
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Table 1. Journal of parallax observations. 


SDSS 

Aref 

-^Vmeas 

N pix 

Epochs 

J0045+1420 

30 

57 

50 

2007.73(4), 2007.82(3), 2008.69(11), 2008.88(16), 2008.97(8), 2009.72(8) 

J0121—0038 

15 

48 

115 

2007.73(8), 2008.05(10), 2008.69(16), 2008.88(17), 2008.97(9), 2009.73(13), 2009.86(18), 
2010.01(10), 2011.75(9), 2011.94(5) 

J0146+1404 

35 

54 

107 

2007.73(8), 2008.05(8), 2008.69(12), 2008.88(18), 2008.97(8), 2009.73(12), 2009.86(16), 
2010.01(12), 2011.75(3), 2011.94(10) 

J0256-0700 

15 

41 

149 

2007.74(33), 2007.81(12), 2008.05(8), 2008.69(13), 2008.88(14), 2009.03(10), 2009.73(8), 
2009.86(12), 2010.02(7), 2011.75(10), 2011.93(22) 

J0301—0044 

25 

58 

102 

2007.73(7), 2007.82(6), 2008.06(8), 2008.69(10), 2008.88(17), 2008.97(8), 2009.73(12), 
2009.86(13), 2010.01(12), 2011.75(9) 

J0309+0025 

17 

47 

126 

2007.74(8), 2007.81(10), 2008.05(8), 2008.69(1), 2008.88(14), 2008.97(7), 2009.72(8), 
2009.86(16), 2010.02(13), 2011.75(16), 2011.93(25) 


{Note. This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 


In most cases our parallaxes were accurate enough that the 
Bayesian adjustments were fairly minor. Furthermore, we 
have four targets in common with the USNO Parallax pro¬ 
gramme and the parallax measurements are in good agree¬ 
ment (Harris, private communication). 

There is only one target in our parallax sample, 
J1547+0523 (NLTT 41210), that does not display signifi¬ 
cant parallax. Th is object was identified as a high proper 
motion target by ILepine fe Sharal (120051) , and included in 
our sample as a WD candidate. We measure relative proper 
motions of —150.5 ± 1.1 and —133.9 ± 1.1 mas yr -1 in RA 
and DEC, respectively. Th ese are consisten t with the proper 
motion measurements by ILepine fe Sharal d2005l ). We also 
measure a parallax of 1.9 ± 1.4 mas, which indicates that 
NLTT 41210 is not a WD. 

2.2 Proper Motion 

In Fig. □ we compare our measured proper motions, as 
listed in Table [2] for the 42 WDs in our sample that also 
have m easured proper mo tions in the SDSS+USNO-B cat¬ 
alogue (M unn c t. al.l 120041) . We expect disagreement at the 
10 mas yr -1 level since our proper motions are relative to 
the particular reference stars used in the reduction. Fig. |T] 
shows that the vast majority of our WDs do indeed fall 
within the range of ± 10 mas yr -1 when compared with the 
SDSS+USNO-B measurements. 

This disagreement arises due to two main factors. First, 
we make no attempt to reduce proper motions to an inertial 
frame. Any systematic trend due to e.g., Galactic rotation 
or solar motion, is still present. Secondly, reference stars 
often have detectable proper motions of their own, so in 
/tra versus fj, dec space they form a cloud of points around 
the origin. Because there are typically only a couple of dozen 
reference stars, the centre of this cloud is statistically uncer¬ 
tain, typically of the order of 5 mas yr -1 . 


2.3 Optical and Infrared Photometry 

We have obtained the available uqriz pho tomet ry from the 
SDSS Data Release 10 (DRIP. lAhn et ahll2014l ) for the 54 
WDs in our sample. These data are listed in columns two 
through six in Table [3] along with their uncertainties. The 


majority of ou r targets al s o hav e near-infrared photometry 
available from lKilic et al.l (|2010l ). and are also listed in Ta¬ 
ble El For the s i x WD s without near-infrared photometry 
from lKilic et al.l (l20ld ) , we adopt the near-infrared photom- 
etry from the UKIDS S Large Area Survey (ULAS) Catalog 
(lLawrence_et_al I[ 2007l) , and the Two Micron All Sky Survey 
(2MASS: lSkrutskie et al.ll2006l) ; see the notes at the bottom 
of Table [3] 


2.4 Optical Spectroscopy 

The majorit y of our targets w e re sel ected from the cool WD 
samples of iKilic et al.l (120061 . l20ld ) , hence they have op¬ 
tical spectroscopy obtained at the McDonald Observatory 
2.7m telescope, Hobby-Eberly Telescope, or the Multiple- 
Mirror Telescope. The ultracool WDs and a few other cool 
WDs have spectroscopy available i n the SDSS or the litera¬ 
ture iLeggett et al.||201 ll: iGiammichele. Bergeron fe Dufouil 


120121 : Tremblay et all 20141 ). There are only eight DA WDs 

in our sample, with the rest of the stars classified as DC 
due to the absence of Ha absorption. This overabundance 
of DC WDs is due to our selection bias for targeting cool 
and ultracool WDs. 


3 THEORETICAL FRAMEWORK 

Our model atmospheres and synthetic spectra are 
derived from the local thermodynamic equilibrium 
(LT E) model atmosphere code o rigina lly described 
in IBereeron. Saumon fe Wesemaell (Il995l ) and refer- 
ences therein, with recent improvements discussed in 
iTremblav fe Bergeronl (120091 ). In particular, we now rely 
on their improved calculations for the Stark broadening 
of hydrogen lines with the inclusion of non-ideal pertur¬ 
bations from protons and electrons - described within the 
occup ation probability formalism of iHummer fe Mihalasl 
(119881 ) directly inside the line profile calculations. Con¬ 
vective energy transport is taken into account following 
the ML2/a = 0.7 prescription of the mixing length theory. 
Non-LTE effects are also included at higher effective 
temperatures but these are irrelevant for the purpose of 
this work. More detai ls regarding our heliu m-atmosphere 
models are provided in iBereeron et al.l (l201ll ) . 
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Table 2. Astrometry of cool WDs. 


SDSS 

RA (J2000) 
(h:m:s) 

Dec. (J2000) 
(d:m:s) 

^abs 

(mas yr -1 ) 

MRA 

(mas yr -1 ) 

mdec 
( mas yr -1 ) 

D 

(pc) 

v ta.n 

(km s -1 ) 

u a 

(km s" 

+ 

V 

(km s“ 

x ) 

W 

(km s - 

- 1 ) 

J0045+1420 

00:45:21.89 

+ 14:20:45.3 

15.9 

± 

1.1 

251 

-52 

eels 

80 

± 

6 

-46 

± 

5 

-37 

± 

4 

-5 

± 

2 

J0121—0038 

01:21:03.00 

-00:38:33.6 

9.6 

± 

1.9 

107 

61 

ns + 

69 

± 

14 

-46 

± 

12 

4 

± 

5 

27 

± 

4 

J0146+1404 

01:46:29.01 

+ 14:04:38.2 

11.0 

± 

1.5 

255 

65 


112 

± 

14 

-82 

± 

13 

-37 

± 

8 

49 

± 

6 

J0256—0700 

02:56:41.62 

-07:00:33.8 

16.4 

± 

1.5 

348 

-193 

eit! 

115 

± 

10 

-17 

± 

3 

-96 

± 

10 

34 

± 

3 

J0301—0044 

03:01:44.09 

-00:44:39.5 

13.4 

± 

1.3 

107 

-538 

+ 6 

192 

± 

16 

83 

± 

7 

-155 

± 

17 

-60 

± 

7 

J0309+0025 

03:09:24.87 

+00:25:25.3 

21.2 

± 

1.6 

-10 

-102 

«!* 

24 

± 

1 

23 

± 

2 

-3 

± 

3 

-4 

± 

2 

J0310—0110 

03:10:49.53 

-01:10:35.3 

7.1 

± 

1.9 

-28 

-77 

164111 

64 

± 

18 

47 

± 

11 

-14 

± 

10 

-25 

± 

10 

J0747+2438N 6 

07:47:21.56 

+24:38:47.7 

18.4 

± 

1.0 

139 

-69 

sell 

40 

± 

2 

29 

± 

2 

-13 

± 

3 

32 

± 

3 

J0747+2438S 6 

07:47:23.50 

+24:38:23.7 

18.4 

± 

1.0 

139 

-69 

eel 1 

40 

± 

2 

29 

± 

2 

-13 

± 

3 

32 

± 

3 

J0753+4230 

07:53:13.28 

+42:30:01.6 

36.2 

± 

1.0 

108 

-388 

27^1 

53 

± 

1 

10 

± 

1 

-40 

± 

2 

10 

± 

1 

J0805+3833 

08:05:57.62 

+38:33:44.1 

47.6 

± 

1.0 

-628 

-552 

2111 

83 

± 

1 

-24 

± 

1 

-30 

± 

1 

-55 

± 

2 

J0817+2822 

08:17:51.52 

+28:22:03.1 

19.6 

± 

1.5 

69 

-197 

5211 

51 

± 

3 

25 

± 

2 

-36 

± 

4 

9 

± 

2 

J0821+3727 

08:21:08.18 

+37:27:38.3 

12.8 

± 

3.2 

167 

-146 

se!?! 

90 

± 

25 

44 

± 

8 

-50 

± 

16 

50 

± 

11 

J0825+5049 

08:25:19.70 

+50:49:20.1 

20.4 

± 

1.3 

-331 

-326 

49! 1 

108 

± 

6 

-48 

± 

4 

-51 

± 

5 

-57 

± 

4 

J0854+3503 

08:54:43.33 

+35:03:52.7 

15.7 

± 

1.5 

-144 

-172 

6411 

68 

± 

6 

-14 

± 

3 

-37 

± 

6 

-32 

± 

4 

J0909+4700 

09:09:14.56 

+47:00:17.5 

16.6 

± 

1.7 

-124 

-172 

eill 

61 

± 

6 

-17 

± 

3 

-36 

± 

6 

-16 

± 

3 

J0942+4437 

09:42:44.96 

+44:37:43.1 

11.7 

± 

1.2 

-144 

-182 

8711° 

96 

± 

9 

-30 

± 

5 

-66 

± 

9 

-23 

± 

4 

J1001+3903 

10:01:03.42 

+39:03:40.5 

11.6 

± 

2.1 

-291 

-175 

67111 

140 

± 

25 

-73 

± 

16 

-75 

± 

16 

-60 

± 

12 

J1107+4855 

11:07:31.38 

+48:55:23.0 

20.9 

± 

1.7 

-730 

-69 

4811 

167 

± 

13 

128 

± 

12 

-56 

± 

6 

-53 

± 

5 

J1115+0033 

11:15:36.97 

+00:33:15.3 

20.2 

± 

2.5 

56 

-246 

sill 

61 

± 

7 

50 

± 

5 

-27 

± 

5 

-14 

± 

3 

J1117+5010 

11:17:08.63 

+50:10:33.9 

21.3 

± 

1.7 

176 

-126 

48l 4 4 

49 

± 

4 

51 

± 

4 

1 

± 

2 

31 

± 

2 

J1158+0004 

11:58:14.52 

+00:04:58.3 

28.9 

± 

1.7 

-25 

183 

sell 

31 

± 

1 

-7 

± 

2 

33 

± 

2 

20 

± 

1 

J1203+0426 

12:03:28.64 

+04:26:53.6 

22.8 

± 

2.1 

-253 

154 

4511 

63 

± 

6 

-50 

± 

6 

13 

± 

2 

10 

± 

1 

J1204+6222 

12:04:39.54 

+62:22:16.4 

18.3 

± 

2.5 

-29 

-154 

58+1° 

43 

± 

6 

16 

± 

2 

-22 

± 

5 

28 

± 

3 

J1212+0440 

12:12:07.01 

+04:40:12.0 

16.0 

± 

2.6 

-280 

-48 

64113 

86 

± 

14 

-54 

± 

11 

-38 

± 

9 

-11 

± 

3 

J1238+3502 

12:38:12.85 

+35:02:49.1 

10.0 

± 

2.0 

-146 

-123 

no+ 26 

100 

± 

20 

-18 

± 

7 

-73 

± 

18 

12 

± 

1 

J1251+4403 

12:51:06.12 

+44:03:03.1 

22.9 

± 

5.3 

22 

-136 

63+ 29 

41 

± 

14 

30 

± 

5 

-8 

± 

5 

15 

± 

2 

J1345+4200 

13:45:32.92 

+42:00:44.2 

27.3 

± 

1.0 

-190 

128 

3611 

40 

± 

1 

-28 

± 

2 

8 

± 

2 

9 

± 

1 

J1349+1155 

13:49:02.33 

+ 11:55:11.8 

35.3 

± 

1.6 

176 

-524 

2811 

74 

± 

3 

69 

± 

3 

-28 

± 

2 

-17 

± 

1 

J1422+0459 

14:22:25.73 

+04:59:39.7 

16.7 

± 

2.1 

-272 

-50 

61 + 

80 

± 

10 

-34 

± 

6 

-49 

± 

8 

28 

± 

3 

J1424+6246 

14:24:29.52 

+62:46:17.1 

21.1 

± 

2.0 

-269 

-42 

4811 

62 

± 

5 

-23 

± 

4 

-31 

± 

5 

33 

± 

3 

J1436+4332 

14:36:42.78 

+43:32:35.7 

37.1 

± 

1.2 

-316 

498 

27+1 

75 

± 

2 

-63 

± 

3 

25 

± 

1 

9 

± 

1 

J1437+4151 

14:37:18.15 

+41:51:51.5 

16.0 

± 

2.2 

-153 

-68 

66+1 2 

52 

± 

7 

-3 

± 

3 

-30 

± 

6 

29 

± 

3 

J1447+5427 

14:47:01.85 

+54:27:44.6 

21.3 

± 

3.5 

-237 

34 

51 + 

58 

± 

10 

-27 

± 

7 

-19 

± 

5 

27 

± 

4 

J1452+4522 

14:52:39.00 

+45:22:38.3 

11.4 

± 

1.4 

-46 

74 

95+ 

39 

± 

5 

-25 

± 

6 

17 

± 

4 

9 

± 

2 

J1458+1146 

14:58:48.52 

+ 11:46:55.9 

17.0 

± 

1.6 

-124 

-96 

soil 

45 

± 

4 

7 

± 

2 

-31 

± 

5 

16 

± 

2 

J1534+4649 

15:34:51.02 

+46:49:49.5 

33.1 

± 

1.8 

-468 

226 

30+ 

75 

± 

4 

-49 

± 

4 

-17 

± 

2 

40 

± 

2 

J1606+2547 

16:06:19.81 

+25:47:02.9 

22.6 

± 

1.7 

-226 

-119 

4511 

54 

± 

4 

13 

± 

2 

-33 

± 

4 

35 

± 

3 

J1615+4449 

16:15:44.67 

+44:49:42.5 

12.1 

± 

3.7 

46 

-231 

89+ 47 
° -20 

99 

± 

37 

100 

± 

27 

-12 

± 

8 

1 

± 

3 

J1632+2426 

16:32:42.23 

+24:26:55.2 

22.9 

± 

1.2 

-12 

-336 

44l 2 2 

70 

± 

3 

65 

± 

3 

-30 

± 

3 

-5 

± 

2 

J1704+3608 

17:04:47.70 

+36:08:47.4 

21.1 

± 

1.7 

183 

-172 

48+ 

57 

± 

4 

52 

± 

4 

20 

± 

2 

-31 

± 

4 

J1722+5752 

17:22:57.78 

+57:52:50.7 

17.8 

± 

1.8 

-37 

395 

5611 

105 

± 

10 

-94 

± 

11 

16 

± 

2 

11 

± 

2 

J1728+2646 

17:28:07.27 

+26:46:19.2 

17.2 

± 

1.8 

-41 

-248 

5911 

70 

± 

7 

67 

± 

6 

-28 

± 

5 

-1 

± 

3 

J2041—0520 

20:41:28.99 

-05:20:27.7 

15.8 

± 

1.6 

-152 

-21 

6511 

47 

± 

5 

41 

± 

4 

3 

± 

3 

41 

± 

4 

J2042+0031 

20:42:59.23 

+00:31:56.6 

16.0 

± 

1.3 

-76 

-241 

6311 

75 

± 

6 

60 

± 

5 

-42 

± 

5 

-8 

± 

3 

J2045—0710 

20:45:57.53 

-07:10:03.5 

12.0 

± 

1.4 

-74 

-125 

sell 2 

59 

± 

7 

49 

± 

5 

-31 

± 

6 

10 

± 

3 

J2118—0737 

21:18:05.21 

-07:37:29.1 

14.5 

± 

2.0 

109 

-127 

72+ 

57 

± 

8 

3 

± 

3 

-23 

± 

6 

-34 

± 

6 

J2147+1127 

21:47:25.17 

+ 11:27:56.1 

18.9 

± 

2.0 

103 

-248 

5411 

69 

± 

7 

26 

± 

3 

-28 

± 

5 

-45 

± 

6 

J2222+1221 

22:22:33.89 

+ 12:21:43.0 

24.4 

± 

1.3 

703 

192 

4ii! 

142 

± 

6 

121 

± 

7 

7 

± 

1 

-44 

± 

3 

J2239+0018A 6 

22:39:54.12 

+00:18:47.3 

12.2 

± 

2.9 

-12 

122 

107+ 

62 

± 

19 

-8 

± 

6 

49 

± 

9 

31 

± 

6 

J2239+0018B b 

22:39:54.07 

+00:18:49.2 

12.2 

± 

2.9 

-12 

122 

1°7+ 

62 

± 

19 

-8 

± 

6 

49 

± 

9 

31 

± 

6 

J2242+0048 

22:42:06.19 

+00:48:22.8 

14.8 

± 

1.7 

133 

-76 

70+ 

51 

± 

5 

-13 

± 

4 

-17 

± 

4 

-24 

± 

4 

J2254+1323 

22:54:08.64 

+ 13:23:57.2 

24.2 

± 

1.6 

326 

-192 

4211 

75 

± 

5 

-24 

± 

3 

-28 

± 

3 

-44 

± 

4 

J2330+0028 

23:30:55.20 

+00:28:52.3 

18.3 

± 

2.3 

137 

104 

5911 

48 

± 

6 

-33 

± 

6 

19 

± 

2 

9 

± 

1 


a Since we do not have any radial velocity measurements for our targets, the U component has been computed assuming i> rac j = 0 km s’ 
k For these two binary systems, a weighted mean was adopted in the determination of their astrometric measurements. 
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Table 3. Optical and near-infrared photometry of cool WDs. 


SDSS 


u 



9 



r 



i 



z 



Y 



J 


H 


K 


J0045+1420 

20.64 

± 

0.08 

19.20 

± 

0.03 

18.45 

± 

0.03 

18.20 

± 

0.03 

18.10 

± 

0.03 


- 


17.24 

± 

0.04 

16.99 

± 

0.04 

16.89 

± 

0.04 

J0121—0038 a 

22.82 

± 

0.28 

20.79 

± 

0.03 

19.74 

± 

0.03 

19.38 

± 

0.03 

19.18 

± 

0.04 

18.47 

± 

0.06 

18.23 

± 

0.08 

18.05 

± 

0.09 

18.10 

± 

0.19 

J0146+1404 

21.37 

± 

0.11 

20.00 

± 

0.03 

19.39 

± 

0.02 

19.27 

± 

0.03 

19.79 

± 

0.11 


- 


19.56 

± 

0.05 

20.07 

± 

0.12 




J0256—0700 

20.74 

± 

0.08 

19.00 

± 

0.02 

18.13 

± 

0.02 

17.79 

± 

0.03 

17.69 

± 

0.03 


- 


16.71 

± 

0.05 

16.62 

± 

0.05 

16.48 

± 

0.06 

J0301—0044 

22.23 

± 

0.34 

20.43 

± 

0.03 

19.38 

± 

0.02 

18.99 

± 

0.02 

18.92 

± 

0.04 


- 


17.96 

± 

0.04 

17.73 

± 

0.04 

17.68 

± 

0.08 

J0309+0025 

19.15 

± 

0.03 

18.19 

± 

0.02 

17.72 

± 

0.02 

17.53 

± 

0.02 

17.50 

± 

0.02 


- 


16.64 

± 

0.04 

16.54 

± 

0.04 

16.87 

± 

0.04 

J0310—0110 

22.71 

± 

0.30 

20.89 

± 

0.04 

20.18 

± 

0.03 

19.91 

± 

0.03 

19.75 

± 

0.08 


- 


18.94 

± 

0.02 

18.73 

± 

0.02 

18.58 

± 

0.02 

J0747+2438N 

21.01 

± 

0.08 

19.29 

± 

0.02 

18.59 

± 

0.02 

18.23 

± 

0.01 

18.14 

± 

0.02 


- 


17.16 

± 

0.04 

16.99 

± 

0.04 

16.85 

± 

0.04 

J0747+2438S 

19.49 

± 

0.03 

18.37 

± 

0.01 

17.91 

± 

0.01 

17.73 

± 

0.01 

17.69 

± 

0.02 


- 


16.78 

± 

0.04 

16.58 

± 

0.04 

16.53 

± 

0.04 

J0753+4230 

19.97 

± 

0.04 

18.09 

± 

0.01 

17.19 

± 

0.01 

16.87 

± 

0.01 

16.75 

± 

0.02 


- 


15.69 

± 

0.04 

15.49 

± 

0.04 

15.47 

± 

0.04 

J0804+2239 

19.73 

± 

0.03 

18.30 

± 

0.02 

17.59 

± 

0.01 

17.39 

± 

0.01 

17.33 

± 

0.02 


- 


16.71 

± 

0.04 

16.92 

± 

0.04 

17.29 

± 

0.06 

J0805+3833 b 

19.00 

± 

0.02 

17.31 

± 

0.01 

16.56 

± 

0.02 

16.27 

± 

0.02 

16.20 

± 

0.02 


- 


15.34 

± 

0.05 

15.19 

± 

0.08 

14.90 

± 

0.09 

J0817+2822 

21.59 

± 

0.16 

19.49 

± 

0.02 

18.61 

± 

0.01 

18.30 

± 

0.01 

18.22 

± 

0.03 


- 


17.33 

± 

0.04 

17.01 

± 

0.04 

16.91 

± 

0.09 

J0821+3727 

20.68 

± 

0.06 

19.14 

± 

0.02 

18.43 

± 

0.01 

18.15 

± 

0.02 

18.04 

± 

0.02 


- 


17.25 

± 

0.04 

17.00 

± 

0.04 

16.85 

± 

0.05 

J0825+5049 

21.09 

± 

0.09 

19.34 

± 

0.02 

18.43 

± 

0.02 

18.09 

± 

0.02 

18.00 

± 

0.03 


- 


17.08 

± 

0.04 

16.83 

± 

0.04 

16.74 

± 

0.04 

J0845+2257 

15.57 

± 

0.01 

15.73 

± 

0.01 

16.08 

± 

0.01 

16.35 

± 

0.02 

16.61 

± 

0.02 


- 


16.24 

± 

0.11 

15.96 

± 

0.00 

16.48 

± 

0.00 

J0854+3503 

23.57 

± 

0.67 

20.53 

± 

0.03 

19.39 

± 

0.02 

19.09 

± 

0.03 

18.95 

± 

0.05 


- 


18.44 

± 

0.04 

18.23 

± 

0.04 

17.98 

± 

0.04 

J0909+4700 

20.64 

± 

0.15 

19.29 

± 

0.03 

18.74 

± 

0.02 

18.50 

± 

0.02 

18.42 

± 

0.05 


- 


18.11 

± 

0.04 

18.62 

± 

0.07 

19.10 

± 

0.10 

J0942+4437 

21.37 

± 

0.09 

19.47 

± 

0.02 

18.58 

± 

0.01 

18.22 

± 

0.02 

18.05 

± 

0.02 


- 


17.15 

± 

0.04 

16.97 

± 

0.04 

16.86 

± 

0.04 

J1001+3903 

21.36 

± 

0.10 

20.05 

± 

0.02 

19.60 

± 

0.02 

20.02 

± 

0.03 

20.61 

± 

0.17 


- 


20.65 

± 

0.06 

21.05 

± 

0.07 




J1107+4855 

21.50 

± 

0.12 

19.49 

± 

0.03 

18.54 

± 

0.02 

18.23 

± 

0.02 

18.11 

± 

0.03 


- 


17.05 

± 

0.05 

16.95 

± 

0.07 

16.86 

± 

0.07 

J1115-|-0033 a 

19.50 

± 

0.04 

17.92 

± 

0.01 

17.22 

± 

0.02 

16.99 

± 

0.01 

16.90 

± 

0.02 


- 


15.78 

± 

0.08 

15.65 

± 

0.18 

15.59 

± 

0.26 

J1117+5010 

21.17 

± 

0.10 

19.34 

± 

0.03 

18.57 

± 

0.03 

18.30 

± 

0.02 

18.16 

± 

0.03 


- 


17.24 

± 

0.04 

17.07 

± 

0.04 

16.97 

± 

0.05 

J1158+0004 

20.86 

± 

0.11 

18.89 

± 

0.04 

17.85 

± 

0.02 

17.54 

± 

0.01 

17.34 

± 

0.03 


- 


16.36 

± 

0.04 

16.31 

± 

0.05 

16.18 

± 

0.05 

J1203+0426 

19.57 

± 

0.03 

18.18 

± 

0.02 

17.50 

± 

0.02 

17.21 

± 

0.01 

17.12 

± 

0.02 


- 


16.39 

± 

0.01 

16.49 

± 

0.02 

16.92 

± 

0.06 

J1204+6222 

20.91 

± 

0.09 

19.25 

± 

0.02 

18.43 

± 

0.02 

18.14 

± 

0.02 

18.06 

± 

0.03 


- 


17.07 

± 

0.04 

16.86 

± 

0.04 

16.80 

± 

0.04 

J1212+0440 

22.07 

± 

0.21 

20.04 

± 

0.03 

19.09 

± 

0.02 

18.79 

± 

0.02 

18.66 

± 

0.04 


- 


17.67 

± 

0.04 

17.50 

± 

0.04 

17.50 

± 

0.05 

J1238+3502 

24.74 

± 

0.81 

21.77 

± 

0.09 

20.31 

± 

0.06 

19.88 

± 

0.05 

20.37 

± 

0.15 


- 


21.19 

± 

0.06 


- 





J1251+4403 

21.46 

± 

0.09 

20.17 

± 

0.03 

20.39 

± 

0.03 

20.72 

± 

0.04 

20.92 

± 

0.17 


- 


21.78 

± 

0.08 


- 





J1345+4200'’ 

19.70 

± 

0.03 

17.85 

± 

0.02 

17.01 

± 

0.01 

16.69 

± 

0.02 

16.54 

± 

0.01 


- 


15.61 

± 

0.06 

15.43 

± 

0.11 

15.00 

± 

0.00 

J1349+1155 

20.55 

± 

0.06 

18.64 

± 

0.02 

17.84 

± 

0.01 

17.42 

± 

0.02 

17.20 

± 

0.02 


- 


16.43 

± 

0.01 

16.29 

± 

0.02 

16.26 

± 

0.02 

J1422+0459 

20.98 

± 

0.10 

19.44 

± 

0.03 

18.58 

± 

0.02 

18.27 

± 

0.02 

18.18 

± 

0.03 


- 


17.15 

± 

0.05 

17.10 

± 

0.08 

17.02 

± 

0.05 

J1424+6246 

20.38 

± 

0.05 

18.83 

± 

0.01 

18.15 

± 

0.03 

17.89 

± 

0.02 

17.71 

± 

0.02 


- 



- 



- 





J1436+4332 

19.83 

± 

0.04 

18.04 

± 

0.02 

17.19 

± 

0.01 

16.85 

± 

0.03 

16.75 

± 

0.03 


- 


15.78 

± 

0.04 

15.62 

± 

0.04 

15.51 

± 

0.04 

J1437+4151 

20.06 

± 

0.04 

19.03 

± 

0.01 

18.45 

± 

0.02 

18.23 

± 

0.03 

18.12 

± 

0.02 


- 


17.43 

± 

0.04 

17.76 

± 

0.05 

18.42 

± 

0.08 

J1447+5427 

21.23 

± 

0.12 

19.46 

± 

0.04 

18.64 

± 

0.02 

18.36 

± 

0.02 

18.25 

± 

0.04 


- 


17.26 

± 

0.07 

17.20 

± 

0.07 

17.07 

± 

0.06 

J1452+4522 

21.55 

± 

0.10 

20.01 

± 

0.02 

19.39 

± 

0.02 

19.31 

± 

0.02 

19.36 

± 

0.06 


- 


18.60 

± 

0.02 

18.43 

± 

0.02 

18.37 

± 

0.02 

J1458+1146 

20.62 

± 

0.08 

18.85 

± 

0.02 

18.02 

± 

0.02 

17.72 

± 

0.02 

17.64 

± 

0.02 


- 


16.63 

± 

0.04 

16.47 

± 

0.05 

16.31 

± 

0.06 

J1534+4649 

20.90 

± 

0.08 

18.76 

± 

0.02 

17.74 

± 

0.02 

17.36 

± 

0.02 

17.19 

± 

0.02 


- 


16.17 

± 

0.04 

16.12 

± 

0.04 

16.04 

± 

0.05 

J1547+0523 

19.96 

± 

0.04 

18.05 

± 

0.01 

17.13 

± 

0.00 

16.75 

± 

0.00 

16.51 

± 

0.01 


- 


15.38 

± 

0.00 

14.95 

± 

0.00 

14.77 

± 

0.01 

J1606+2547 

20.99 

± 

0.08 

19.24 

± 

0.02 

18.45 

± 

0.02 

18.17 

± 

0.02 

18.07 

± 

0.04 


- 


17.07 

± 

0.04 

17.09 

± 

0.06 

16.84 

± 

0.06 

J1615+4449 

21.18 

± 

0.10 

19.59 

± 

0.02 

18.84 

± 

0.02 

18.57 

± 

0.02 

18.52 

± 

0.04 


- 


17.44 

± 

0.04 

17.24 

± 

0.05 

17.26 

± 

0.07 

J1632+2426 

21.47 

± 

0.10 

19.60 

± 

0.02 

18.73 

± 

0.02 

18.49 

± 

0.02 

18.40 

± 

0.03 


- 


17.67 

± 

0.02 

18.10 

± 

0.02 

18.04 

± 

0.02 

J1704+3608 

20.50 

± 

0.05 

18.72 

± 

0.01 

17.94 

± 

0.01 

17.66 

± 

0.01 

17.55 

± 

0.02 


- 


16.62 

± 

0.04 

16.34 

± 

0.04 

16.32 

± 

0.06 

J1722+5752 

20.39 

± 

0.06 

19.28 

± 

0.02 

18.79 

± 

0.02 

18.56 

± 

0.03 

18.50 

± 

0.03 


- 


17.74 

± 

0.04 

17.84 

± 

0.05 

18.75 

± 

0.12 

J1728+2646 

19.18 

± 

0.03 

18.14 

± 

0.02 

17.68 

± 

0.01 

17.51 

± 

0.01 

17.46 

± 

0.02 


- 



- 



- 





J2041—0520 

20.95 

± 

0.08 

19.27 

± 

0.01 

18.51 

± 

0.01 

18.24 

± 

0.01 

18.14 

± 

0.03 


- 


17.25 

± 

0.04 

16.97 

± 

0.04 

16.93 

± 

0.04 

J2042+0031 

21.67 

± 

0.14 

19.95 

± 

0.02 

19.05 

± 

0.01 

18.73 

± 

0.01 

18.61 

± 

0.03 


- 


17.65 

± 

0.04 

17.45 

± 

0.04 

17.36 

± 

0.05 

J2045—0710 

21.03 

± 

0.09 

19.33 

± 

0.02 

18.60 

± 

0.01 

18.33 

± 

0.01 

18.15 

± 

0.03 


- 


17.32 

± 

0.04 

17.10 

± 

0.04 

17.03 

± 

0.04 

J2118—0737 

23.38 

± 

0.95 

20.70 

± 

0.03 

19.48 

± 

0.03 

19.01 

± 

0.02 

18.76 

± 

0.04 


- 


17.90 

± 

0.04 

17.82 

± 

0.04 

17.81 

± 

0.05 

J2147+1127 

20.83 

± 

0.09 

19.19 

± 

0.02 

18.43 

± 

0.02 

18.13 

± 

0.02 

18.01 

± 

0.03 


- 


17.14 

± 

0.04 

16.84 

± 

0.04 

16.79 

± 

0.04 

J2222+1221 

21.74 

± 

0.11 

19.48 

± 

0.02 

18.37 

± 

0.02 

17.88 

± 

0.02 

17.67 

± 

0.02 


- 



- 



- 





J2239+0018A a 

21.27 

± 

0.08 

20.14 

± 

0.04 

19.59 

± 

0.03 

19.48 

± 

0.05 

20.28 

± 

0.16 

19.57 

± 

0.10 

19.69 

± 

0.19 


- 





J2239+0018B a 

23.13 

± 

0.36 

20.79 

± 

0.04 

19.88 

± 

0.03 

19.49 

± 

0.03 

19.24 

± 

0.06 

18.66 

± 

0.05 

18.34 

± 

0.06 

17.98 

± 

0.10 

18.48 

± 

0.27 

J2242+0048 

22.11 

± 

0.22 

19.63 

± 

0.02 

18.65 

± 

0.02 

18.28 

± 

0.01 

18.16 

± 

0.03 


- 


18.06 

± 

0.04 

18.72 

± 

0.07 

19.16 

± 

0.10 

J2254+1323 

21.57 

± 

0.17 

19.51 

± 

0.02 

18.49 

± 

0.02 

18.14 

± 

0.01 

18.00 

± 

0.02 


- 


17.04 

± 

0.04 

16.88 

± 

0.04 

16.85 

± 

0.04 

J2330+0028 

21.85 

± 

0.24 

19.88 

± 

0.02 

18.95 

± 

0.02 

18.66 

± 

0.02 

18.53 

± 

0.04 


- 


17.63 

± 

0.04 

17.36 

± 

0.04 

17.32 

± 

0.04 


a IR photometry from UKIDSS 
b IR photometry from 2MASS 


Our model grid covers a range of effective temperature 
between T e g = 1500 and 45,000 K in steps of 500 K for T e g 
< 15,000 K, 1000 K up to T eS = 18,000 K, 2000 K up to T eS 
— 30,000 K and by steps of 5000 K above. The logg ranges 
from 6.5 to 9.5 by steps of 0.5 dex, with additional models 
at log g = 7.75 and 8.25. We also calculated mixed hydrogen 
and helium atmosphere models with log (He/H) = —2.0 to 
5.0, in steps of 1.0 dex. 

Since the photometric technique described below relies 
heavily on the flux at the u and g bandpasses, we now in- 
clude in our models the op acity from the red wing of Lya 
dKowalski fe Saumonl 1200^ 1. which significantly affects the 
flux in the ultraviolet. 


4 PHOTOMETRIC ANALYSIS 

4.1 General Procedure 

Atmospheric parameters, T e fj and logg, and chemi¬ 
cal compositions of cool WDs can be measured ac- 
curately using the photometric technique developed by 
iBergeron, Ruiz fe Leggett] (Il997l l. We first convert optical 
and infrared photometric measurements into observed fluxes 
and compare the resulting energy distributions with those 
predicted from our model atmosphere calculations. To ac¬ 
complish this task, we first transform every magnitude m 
into an average flux f™. Since ugriz photometry is defined 
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Figure 2. Location of the WDs in our sample in a M r versus r — i (left) and r — z (right) colour-magnitude diagram. The black dots 
correspond to the WDs with H-rich atmospheres while the white dots represent the He-rich WDs. The red dots represent the 15 cool 
and ultracool WDs with mixed atmospheres. The solid and dashed black lines represent pure H and pure He tracks, respectively, for 
masses M = 0.3, 0.6, and 0.9 Mq, from right to left. The solid blue lines represent the predictions from mixed model atmospheres for 
M = 0.2 Mq. The mixed atmosphere model tracks are labelled with their He abundance log (He/H) = —2 and 2. LHS 3250 and J1102 
are shown as a green triangle, and a blue square, respectively. 


on the AB magnitude system, we first calculate f™ using 
the equation 


m = —2.5 log f™ — 48.60 (1) 

and then f™ is converted to f™ following f ™ = f™c/\ 2 , 
where A is the central wavelength of the given filter. For the 
near-infrared photometry, we obtain using the equation 


m =-2.5 log f™ + Cm, (2) 

where c m is a constant to be determined for each filter, as 
described below. In general, 


are determined usin g the improved calibration fluxes from 
iHolberg fc Bergeron ll20 06ll. defined w ith the HST absolute 
flux scale of Vega ( Bohlin fe Gillila nd 2004), and appropri¬ 
ate magnitudes on a given system. 

For each star in Table [3] a minimum set of five average 
fluxes f™ is obtained, which can be compared with model 
fluxes. Since the observed fluxes correspond to averages over 
given bandpasses, the monochromatic fluxes from the model 
atmospheres need to be converted into average fluxes, H 
by substituting f\ in Equation © for the monochromatic 
Eddington flux, H\. We can then relate the average observed 
fluxes f™ and the average model fluxes H™ - which depend 
on Teff, log grand chemical composition - by the equation 


f 0 °° fxSm(\)\d\ 

Ix / 0 °° S m (A)AdA W 

where S m ( A) is the transmission function of the corre¬ 
sponding bandpass, f\ is the monochromatic flux from 
the star received at Earth. For the ugriz photometry, a 
slightly different definition o f the above Equation © is re- 
quired (see Equation (3) of IHolberg fc Bergeronl fiood . for 
instance). The trans mission functions for the u c/riz sys¬ 
tem are described in IHolberg fe Bergeronl (l2006l l and ref¬ 
erences therein. The transmission functions for the JHK or 
JH Ks filters on the MKO photometr ic system are taken 
from lTokunaga, Simons &: Vaccal l2002h . 

The constants c m in Equation © for each passband 


f? = 4 n(R/D) 2 H? (4) 

where R/D defines the ratio of the radius of the star to 
its distance from Earth. We then minimize the \ 2 value 
defined in terms of the difference between observed and 
model fluxes over all bandpasses, properly weighted by the 
photometric uncertainties. Our minimization procedure re¬ 
lies on the non-lin ear least-s quares method of Levenberg- 
Marquardt (IPress et al.lll986l l. which is based on a steepest 
decent method. Only T e g and the solid angle n(R/D) 2 are 
considered free parameters, while the uncertainties of both 
parameters are obtained directly from the covariance matrix 
of the fit. 

For stars with known trigonometric parallax mea¬ 
surements, we first assume a value of log<? = 8.0 and 
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determine the effective temperature and the solid an¬ 
gle, which combined with the distance D obtained from 
the trigonometric parallax measurement, yields directly 
the radius of the star R. The radius is then converted 
into mass using evolutionary models sim il ar to those de¬ 
scribed in iFontaine. Brassard fe Bergeronl (l200lh but with 
CO cores, g(He) = logMHe/M* = 10 -2 and q(H) = 10 -4 , 
which are representative of hydrogen-atmosphere WDs, and 
q( He) = ICC 2 and g(H) = 10~ 10 , which are representa¬ 
tive of helium-atmosphere WDs. After the first iteration, 
if M < 0.406 Mb, we switch to the ev olutionary models of 
lAlthaus, Serenelli fe Benvenutol (l200ll l , appropriate for low- 
mass He-core WDs. In general, the log g value obtained from 
the inferred mass and radius (g = GM/R 2 ) will be different 
from our initial guess of loggr = 8.0, and the fitting proce¬ 
dure is thus repeated until an internal consistency in logf/ 
is reached. 


4.2 Results 


Fig. [5] presents the colour-magnitude diagram for our par¬ 
allax sample along with the evolutionary tracks for 0.3- 
0.9 Mq pure H, pure He, and 0.2 Mg mixed H/He atmo¬ 
sphere models. Note that all the evolutionary tracks plot¬ 
ted in Fig. [2] represent the evolution of CO-core WDs. 
Two other ultracool WDs with parallax measurements and 
SPSS photometry, LHS 32 50 an d J1102 ([Harris et al.l 19991 : 
Ber geron, Leggett fc RuI3l200ll ; lHall et ah 20081 : Kilic et ahl 
2012|), are also included for comparison. 

Interestingly, the majority of the targets in our sample 
fall above the evolutionary tracks for 0.6 Mq WDs, indi¬ 
cating that they are low-mass objects. Some of these WDs 
are even brighter than the 0.3 Mq WD sequence, implying 
masses as low as « 0.2 Mq. A significant fraction of the 
stars in our sample are IR-faint WDs that suffer from CIA 
from molecular hydrogen. The CIA affects the redder opti¬ 
cal bands and the infrared. Hence, most of these IR-faint 
objects lie to the left of the pure H and pure He model se¬ 
quences. Note that our sample was selected to include as 
many IR-faint WDs as possible. Therefore, these are over¬ 
represented in this figure. It is clear from this figure that the 
colour-magnitude distribution of our sample is well matched 
by WD models with masses « 0.2-0.9 Mq with a variety of 
compositions, including pure H, pure He and mixed H/He 
atmospheres. Below we discuss the DA, DC and ultracool 
WD samples separately. 


4.2.1 DA WDs 

Fig. [3] displays the best-fitting pure-hydrogen models to 
the SEDs of the eight WDs classified as DA. Both the ob¬ 
served SEDs and the Ho line profiles are reproduced fairly 
well by our pure H models. Given our parallax measure¬ 
ments, the best-fitting radii for these eight targets range 
from 0.011 to 0.022 Rq (R > R0), indicating that they 
are relatively low-mass WDs. In fact, half of these WDs 
have masses below 0.45 Mq, and therefore are likely He-core 
WDs. The majority of low- mass WDs are in short-period 
(P < 1 d) binary systems llMarsh. Dhillon Duck|ll995l : 
iBrown et all l201l|j . Therefore, J0045+1420, J0821+3727, 
.11115+0033, and J1728+2646 are likely unresolved bi¬ 
nary WDs. Table 4 provides WD cooling age estimates 


for these DA WDs, as well as the rest of our paral¬ 
lax sample. For M < 0.45 Mq WDs, we provide cool¬ 
ing ages for both CO and He core composition based on 
the ev olutiona ry tracks of Fontain^^rassard^&Bergeron) 
(120011 1 and lAlthaus. Serenelli fe Benvenutol ( 200 il l. respec¬ 
tively. Regardless of the core composition, these eight DA 
WDs have cooling ages of less than 8 Gyr. 

It is necessary to note an important caveat regarding 
the four potential binaries listed above. If they are indeed 
unresolved binaries, then the WDs in these systems will be 
more massive than implied by our fits assuming a single star. 
Hence, their actual cooling ages will be larger for a given T e g. 
Our estimates for the cooling ages of these potential binaries 
should therefore be regarded, at best, as lower limits. 


4.2.2 DC WDs 


Fig. H shows our model fits to the SEDs of the 31 DC WDs 
that are best explained by pure H or pure He atmosphere 
models. In all cases, the optical spectra are featureless near 
the Ha region. Hence, the choice of a pure H or pure He 
composition is based solely on the fits to the optical and 
infrared photometry. In most cases, the atmospheric param¬ 
eters from both the pure H and pure He solution agree 
within the uncertainties. Our model fits indicate that all 
of these WDs have T e « < 5000 K. The ratio of the H to 
He atmosphere WDs is 13/18. However, all DC WDs with 
temperatures below T e g = 4530 K are best explained by H- 
rich atmosphere models (see also Kowalski fe Saumonll2006l : 


iGiammichele, Bergeron fc Dufourll2012l')~ 

Just like the DA sample discussed above, about half of 
the DCs in this sample are low-mass objects. The two coolest 
stars, J2118—0737 and J2222+1221, have T e g = 3920 ± 60 
and 4010 ± 80 K, and M = 0.31 ± 0.09 Mq and 
0.37 ± 0.03 Mq, respectively. Assuming He-cores, these tem¬ 
peratures correspond to cooling ages of 7.7 and 9.4 Gyr, re¬ 
spectively. If these are short-period, unresolved binary sys¬ 
tems, then the companions would be fainter and more mas¬ 
sive WDs. Due to the unknown prior history of such binary 
systems and without an estimate on their initial masses, 
their total ages, including the main-sequence + WD cooling 
ages, cannot be reliably calculated. 


4-2.3 DC WDs with Mixed H/He Atmospheres 

iGates et al.l (120041 1 . lHarris et al.l (l2008l l and iKilic et~ahl 
( 2010l l have identified several IR-faint WDs that were orig¬ 
inally thought to be ultracool WDs with T e g < 4000 K. It 
turns out that some of these IR-faint WDs are relatively 
warm. There are nine IR-faint, DC WDs in our sample that 
are best-fitted with T e g > 4500 K mixed H/He atmospheres 
models. The main opacity source in these mixed models is 
the H 2 -He CIA in the infrared. Since cool He-rich WDs have 
lower opacities and higher atmospheric pressures, the CIA 
becomes effective at high er temperatures (T e g > 4000 K, 
iBergeron fe Leggettll2002l l. 

Fig. [5] shows the SEDs for these nine DC WDs with 
mixed composition. The mixed models with log(He/H) = 
—0.4 to 2.3 fit the observed SEDs (over the 0.3-2.2 /im re¬ 
gion) fairly well. The best-fitting parameters for some of 
these sta rs are mark edly different than the parameters pre¬ 
sented in lKilic et al.l J201ol l. However, the analysis presented 
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Figure 3. Fits to the observed energy distributions with pure hydrogen models (filled circles) and with pure helium models (open 
circles) for the eight WDs that exhibit, or potentially exhibit, absorption at Ha. Adopted atmospheric parameters are emphasized in red. 
Here and in the following figures, the photometric observations are represented by error bars while the filled and open circles represent 
the model fluxes for the pure H and pure He solution, respectively. In the right-hand panels we show the observed normalized spectra 
together with the synthetic line profiles calculated with the atmospheric parameters corresponding to the pure hydrogen solutions. 


in this paper is superior to earlier work since we now in¬ 
clude all available photometry in our analysis (including the 
u-band data) and we also have parallax measurements avail¬ 
able. J1632+2426 is the most-massive and the oldest WD (in 
terms of the WD cooling age) in this sample, with a mass 
of 0.82 ± 0.04 Mq and a cooling age of 7.7 Gyr. 


4-2.4 Ultracool WDs 

We originally selected 12 ultracool WD candidates 
for follow-up parall ax observat io ns: J 0854+3503 and 
J1001+3903 from iGates et al.l (120041) ; J0121-0038, 

J0301—0044, J2239+0018 and J2242+0048 from 

IVidrih et al.l J2007l) : J0146+1404, J0310-0110, 

J123 8+3502, J1251+440 3, J1452+4522 and J1632+2426 
from [Harris et all d2008l ). Our detailed model atmosphere 
analysis using parallax data shows that only half of these 
stars are actually ultracool WDs with T e g < 4000 K. The 
rest of the ultracool candidates are best explained by pure 
H/He or mixed atmosphere models with T e g > 4000 K. 

Fig. [S] shows the SEDs and our model fits to the six ul¬ 
tracool WDs in our sample. The best-fitting parameters are 
given in each panel and at the end of Table 4. Note that prior 
to this work, there were only three ultracool WDs with par¬ 
allax observations available. Hence, the ultracool WD sam¬ 
ple presented here is a significant addition to this sample. 
The six ultracool WDs presented here are best explained 
by mixed H/He atmospheres with T e fj = 2710-3760 K and 
log(He/H) = 0.65-2.96. Interestingly, all six of these ultra¬ 


cool WDs are too bright for average mass WDs. Instead, the 
observed parallaxes require relatively large radii (R = 0.015- 
0.023 Rq) and low masses (M = 0.17-0.39 Mq). Assuming 
He-cores, the WD cooling ages range from 4.5 to 9.7 Gyr. 
They are located within 63-110 pc of the Sun and they dis¬ 
play tangential velocities of 40-140 km s -1 . Hence, these 
ultracool WDs likely belong to the Galactic disc. 


5 DISCUSSION 
5.1 Nearby WDs 

The local WD population is complete to within 13 pc, 
and there remains a significant numb er of WDs to be dis¬ 
covered in the solar neighbourhood llHolberg et al.l 120081 : 
iGiammichele, Bergeron fe DufouiH2012l ). Through our par¬ 
allax observations, here we have uncovered WDs w ith dis¬ 
tances ranging from 21 to «100 pc. lSion et al.l (120141 ) present 
224 WDs within 25 pc of the Sun. With a distance of 
21 ± 1 pc, J0805+3833 (WD 0802+387) is a new addition 
to this sample. There are also four other WDs, J0753+4230, 
J1349+1155, J1436+4332 and J1534+4649, with distances 
^ 30 pc. Since parallax observations on individual targets is 
time consuming, significant progress on creating a complete 
sample of WDs in the solar neighbourhood has to wait un- 
til astrometric data fr om large scale surveys such as GAIA 
(IPerrvman et al. ]|200l|) a nd the Large Synoptic Survey Tele- 
scope (LSST: llvezic et al.ll2008l ) become available. 
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Figure 4. Fits to the observed energy distributions with pure hydrogen models (filled circles) and with pure helium models (open 
circles) for the 31 DC WDs. All objects have featureless spectra near the Ho; region, and the SEDs are best explained with pure model 
atmospheres. Adopted atmospheric parameters are emphasized in red. 
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Figure 4. - continued 
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Figure 5. Fits to the SEDs of the nine IR-faint, DC WDs in our sample, excluding the ultracool WDs. All objects have featureless 
spectra near the Ha region, and the SEDs are best explained with mixed model atmospheres. Note that the measured abundances are 
quoted relative to the dominant atmospheric constituent. 
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Table 4. Properties of cool WDs. 


SDSS 

T eff 

logs 

M 

R 

Comp. 

Mg 

T cool,CO 

7"cool,He 


(K) 

(cm s -2 ) 

(Mq) 

(R 0 ) 

(log He/H) 

(mag) 

(Gyr) 

(Gyr) 

DA 


J0045+1420 

5090 

± 

60 

7.73 ± 0.15 

© 

CO 

H- 

0.07 

0.015 ± 0.001 

H 

15.10 ± 

0.18 

9 7 +0.9 
^ - 0.6 

- 

J0747+2438S 

5590 

± 

40 

7.92 =1= 0.09 

0.54 ± 

0.05 

0.013 ± 0.001 

H 

14.67 ± 

0.12 

2 5 + 0 ' 5 
—0.3 

- 

J0821+3727 

5050 

± 

50 

7.27 ± 

0.50 

0.27 ± 

0.15 

0.020 ± 0.006 

H 

14.47 ± 

0.59 

1 6 +1 -° 

—1.3 

3 1 + 24 
cil -l.l 

J1115+0033 

4910 

± 

40 

7.05 ± 0.27 

0.21 

± 

0.07 

0.022 ± 0.003 

H 

14.38 ± 

0.28 

1 4 + 0 - 4 

1 -0.7 

9 q+ 0.6 
^• y - 0.6 

J1424+6246 

4970 

± 

60 

7.88 ± 0.16 

0.51 

± 

0.09 

0.014 ± 0.001 

H 

15.42 

± 

0.20 

4 4+1-7 

- 

J1606+2547 

4860 

± 

60 

8.18 ± 0.12 

0.69 ± 

0.08 

0.011 ± 0.001 

H 

15.97 ± 

0.17 

7 q+ 0-6 

- 

J1728+2646 

5600 

± 

50 

7.65 ± 0.17 

0.42 

± 

0.07 

0.016 ± 0.002 

H 

14.29 

± 

0.24 

1 6 +0 ' 4 
xo — 0.2 

o c+ 0.8 
°-°- 0.6 

J2147+1127 

4920 

± 

60 

7.90 ± 0.18 

0.52 ± 

0.10 

0.013 

± 0.002 

H 

15.53 

± 

0.24 

4 9 + 1 - 8 
*‘ y — 1.8 

- 

DC 

J0121 —0038 

4560 

± 

50 

7.31 

± 

0.39 

0.28 

± 

0.13 

0.019 

± 

0.004 

He 

15.43 

± 

0.46 

2 2"*” 1 ' 1 
zz -1.2 

4 6+ 3 ' 2 
— 1.6 

J0256—0700 

4780 

± 

40 

7.39 

± 

0.15 

0.30 

± 

0.05 

0.018 

± 

0.002 

He 

15.07 

± 

0.20 

9 n+0.3 
ZZ -0.3 

4 6"*” 1 ' 4 
*‘° —0.9 

J0301—0044 

4530 

± 

50 

7.79 

± 

0.16 

0.45 

± 

0.08 

0.014 

± 

0.001 

He 

16.08 

± 

0.19 

4 8 +1 ' 6 
° — 1.2 

- 

J0310—0110 

4990 

± 

50 

7.48 

± 

0.47 

0.34 

± 

0.17 

0.018 

± 

0.005 

He 

14.82 

± 

0.59 

2 1+ 21 

4 7 + 3 ’ 2 
-2.0 

J0747+2438N 

4780 

± 

40 

7.83 

± 

0.10 

0.47 

± 

0.05 

0.014 

± 

0.001 

He 

15.59 

± 

0.12 

4 5 +11 
*'° —0.9 

- 

J0753+4230 

4580 

± 

40 

7.80 

± 

0.05 

0.46 

± 

0.03 

0.014 

± 

0.001 

H 

15.89 

± 

0.06 

C 1 + 0.6 
D, 1 - 0.6 

- 

J0805+3833 

4890 

± 

40 

8.02 

± 

0.03 

0.59 

± 

0.02 

0.012 

± 

0.001 

H 

15.70 

± 

0.05 


- 

J0817+2822 

4660 

± 

50 

7.90 

± 

0.13 

0.51 

± 

0.08 

0.013 

± 

0.001 

He 

15.91 

± 

0.17 

5-6±i | 

- 

J0825+5049 

4660 

± 

40 

7.86 

± 

0.11 

0.49 

± 

0.06 

0.014 

± 

0.001 

He 

15.89 

± 

0.14 

c 9 + 1 ' 3 

1.0 

- 

J0942+4437 

4770 

± 

40 

7.10 

± 

0.21 

0.22 

± 

0.06 

0.022 

± 

0.002 

He 

14.77 

± 

0.23 


o q + 0.6 

**-0.5 

J1107+4855 

4510 

± 

70 

7.88 

± 

0.14 

0.50 

± 

0.08 

0.014 

± 

0.001 

H 

16.08 

± 

0.18 


- 

J1117+5010 

4870 

± 

60 

8.17 

± 

0.12 

0.69 

± 

0.08 

0.011 

± 

0.001 

H 

15.93 

± 

0.18 

7 o+ 0.8 

7 - 8 -l.l 

- 

J1158+0004 

4320 

± 

70 

7.72 

± 

0.11 

0.42 

± 

0.05 

0.015 

± 

0.001 

H 

16.16 

± 

0.13 

K 9+1 ■ 1 
O- ^- 1.0 

- 

J1204+6222 

4800 

± 

40 

7.73 

± 

0.23 

0.44 

± 

0.11 

0.015 

± 

0.003 

He 

15.43 

± 

0.31 

7 2 + 4 ' 3 
‘ Z -2.0 

- 

J1212+0440 

4530 

± 

80 

7.85 

± 

0.29 

0.49 

± 

0.16 

0.014 

± 

0.003 

H 

16.01 

± 

0.37 

K Q + 2 - 5 
— 2.8 

- 

J1345+4200 

4690 

± 

20 

7.25 

± 

0.06 

0.26 

± 

0.02 

0.020 

± 

0.001 

He 

15.02 

± 

0.08 

2 0+ 01 

4 0 +O - 3 
*' u —0.3 

J1422+0459 

4800 

± 

50 

7.71 

± 

0.24 

0.41 

± 

0.11 

0.015 

± 

0.002 

He 

15.51 

± 

0.29 

o r + 2.2 

^■°-i.o 

- 

J1436+4332 

4630 

± 

30 

7.88 

± 

0.05 

0.50 

± 

0.03 

0.013 

± 

0.001 

He 

15.88 

± 

0.07 

5 ' 4 io ,6 

- 

J1447+5427 

4750 

± 

60 

8.04 

± 

0.30 

0.59 

± 

0.18 

0.012 

± 

0.003 

He 

15.92 

± 

0.40 

6-9±3.1 

- 

J1458+1146 

4810 

± 

40 

7.35 

± 

0.15 

0.29 

± 

0.05 

0.019 

± 

0.002 

He 

14.96 

± 

0.20 

9 -,+0.3 

1 — 0.3 

4 2 + 1 ' 3 
*'^ — 0.8 

J1534+4649 

4270 

± 

60 

7.82 

± 

0.09 

0.47 

± 

0.05 

0.014 

± 

0.001 

H 

16.35 

± 

0.12 

6 4 + 1 - 0 

- 

J1615+4449 

4930 

± 

60 

7.41 

± 

0.57 

0.31 

± 

0.18 

0.018 

± 

0.007 

H 

14.84 

± 

0.74 

2.0t\i 

4 3 + 3 ’ 3 
—2.3 

J1704+3608 

4760 

± 

40 

7.64 

± 

0.13 

0.40 

± 

0.06 

0.016 

± 

0.001 

He 

15.31 

± 

0.18 

9 q + 0.9 
2 - y -0.4 

6 . 8 iU 

J2041—0520 

4870 

± 

40 

7.67 

± 

0.17 

0.42 

± 

0.07 

0.016 

± 

0.002 

He 

15.20 

± 

0.23 

9 o + l • 1 
2 - 8 - 0.6 

6 4 + 1 - 5 
° -1.3 

J2042+0031 

4660 

± 

50 

7.91 

± 

0.15 

0.52 

± 

0.09 

0.013 

± 

0.001 

He 

15.95 

± 

0.19 

c 7+1.5 
b -‘ - 1.6 

- 

J2045—0710 

4950 

± 

40 

7.30 

± 

0.21 

0.28 

± 

0.06 

0.019 

± 

0.003 

He 

14.66 

± 

0.26 

1 s +°- 4 
1 , 8 -0.4 

3 5 + 1,4 
°-°- 0.6 

J2118—0737 

3920 

± 

60 

7.42 

± 

0.25 

0.31 

± 

0.09 

0.018 

± 

0.003 

H 

16.41 

± 

0.32 

3 4 + 1- 8 
o -*- 0.8 

7 7+3.6 

'• ‘ -2.3 

J2222+1221 

4010 

± 

80 

7.56 

± 

0.08 

0.37 

± 

0.03 

0.017 

± 

0.001 

H 

16.42 

± 

0.11 

4 9 + 0 - 9 
* - 0.6 


J2239+0018B 

4420 

± 

90 

7.46 

± 

0.55 

0.33 

± 

0.19 

0.018 

± 

0.006 

H 

15.64 

± 

0.66 

9 q+3-0 
zy - 2.2 


J2254+1323 

4350 

± 

70 

7.95 

± 

0.12 

0.54 

± 

0.07 

0.013 

± 

0.001 

H 

16.39 

± 

0.16 


- 

J2330+0028 

4650 

± 

80 

7.98 

± 

0.22 

0.57 

± 

0.13 

0.013 

± 

0.002 

H 

16.03 

± 

0.29 

6 o+l - 7 
—2.7 

- 

Mixed 

J0309+0025 

5610 

± 

70 

7.98 

± 

0.13 

0.56 

± 

0.08 

0.013 

± 

0.001 

1.30 

14.74 

± 

0.18 

3 4+1-1 
° - 0.8 

- 

J0854+3503 

4550 

± 

80 

8.22 

± 

0.15 

0.71 

± 

0.10 

0.011 

± 

0.001 

-0.23 

16.50 

± 

0.22 

7 6+ 0 - 2 

'•°-o.i 

- 

J0909+4700 

4510 

± 

110 

7.79 

± 

0.19 

0.45 

± 

0.10 

0.014 

± 

0.002 

2.31 

15.36 

± 

0.23 

4 8+ 2 ' 0 
*’° —1.4 

- 

J1203+0426 

5010 

± 

40 

7.65 

± 

0.16 

0.41 

± 

0.07 

0.016 

± 

0.002 

1.08 

14.91 

± 

0.22 

5 6+ 1 ' 0 
°-°- 1.0 

- 

J1349+1155 

4710 

± 

40 

8.30 

± 

0.06 

0.76 

± 

0.04 

0.010 

± 

0.001 

-0.40 

16.37 

± 

0.10 

7 5 + 01 
'• b - 0.1 

- 

J1437+4151 

5000 

± 

70 

7.75 

± 

0.24 

0.45 

± 

0.11 

0.015 

± 

0.003 

1.97 

14.93 

± 

0.32 


- 

J1452+4522 

5580 

± 

70 

8.25 

± 

0.20 

0.73 

± 

0.13 

0.011 

± 

0.002 

0.71 

15.12 

± 

0.30 

5 - 5 -l.t 

- 

J1632+2426 

4650 

± 

70 

8.39 

± 

0.06 

0.82 

± 

0.04 

0.010 

± 

0.001 

0.55 

16.38 

± 

0.10 

77 + 0 .l 
'•' - 0.1 

- 

J1722+5752 

5230 

± 

80 

8.34 

± 

0.14 

0.79 

± 

0.09 

0.010 

± 

0.001 

1.49 

15.54 

± 

0.21 

6 7 + 0 ' 1 
°‘ -0.5 

- 


Ultracool 


J0146+1404 

3600 

± 

140 

7.59 ± 0.22 

0.38 ± 0.09 

0.016 ± 0.003 

2.82 

15.23 ± 0.29 

3 7 + 1 - 7 
^•'- 1.0 

®-a±i:S 

J1001+3903 

2710 

± 

150 

7.06 ± 0.24 

0.20 ± 0.09 

0.022 ± 0.004 

2.93 

15.35 ± 0.39 

9 q+2.3 
2 - y -0.7 

7.6ti- 2 7 

J1238+3502 

2900 

± 210 

7.13 ± 0.42 

0.22 ± 0.11 

0.021 ± 0.004 

0.65 

16.77 ± 0.44 

3 6+ 3 ' 2 

q 7 + 6-8 
y - '-5.0 

J1251+4403 

2750 

± 

180 

7.67 ± 0.69 

0.39 ± 0.28 

0.015 ± 0.006 

2.91 

16.17 ± 0.73 

9 1 + 1 ' 4 
y ‘ 1 — 4.9 

9-i±i:S 

J2239+0018A 

3510 

± 220 

6.95 ± 0.34 

0.20 ± 0.09 

0.025 ± 0.008 

2.88 

14.99 ± 0.66 

2 0 +1 - 8 
2 U -0.9 

™±i:l 

J2242+0048 

3770 

± 

90 

7.00 ± 0.13 

0.20 ± 0.03 

0.023 ± 0.003 

0.89 

15.41 ± 0.25 

i q+ 0.8 
J " y — 1.2 
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Figure 6. Fits to the SEDs of the six ultracool DC WDs in our sample. All objects have featureless spectra near the Ho; region, and 
the SEDs are best explained with mixed model atmospheres. Note that the measured abundances are quoted relative to the dominant 
atmospheric constituent. 
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Figure 7. Similar to Fig.[6]but for the ultracool WD LHS 3250. Figure 8 . Similar to Fig.[4]but for the ultracool WD J1102+4113 

where we show the pure H solution only. 


5.2 The Nature of Ultracool WDs 


The colour-magnitude diagram presented in Fig. [2] shows 
that the faintest WDs in our sample have M r ^ 16.4 mag. 
This magnitude limit corresponds to WD cooling ages of 
11 Gyr for 0.6 Mq CO-core and pure H atmosphere WDs. 
There are six ultracool WDs in our sample with best-fitting 
temperatures of < 4000 K. Constraining the nature, includ¬ 
ing the total ages of th ese stars, has been p roblematic. 

lOppenheimer et all 


Previously, _ lOppenhei mer et all _ (|2001), 

iBergeron, Leggett fe Ruiz ( 20011 1. and Kilic et all (1 201 21 1 


presented detailed model atmosphere analysis of three 
ultracool WDs with parallax measurements. WD 0346+246 
is a 3650 K, logt/ = 8.3 mixed atmosphere WD with 
log(He/H) = 0.43, whereas J1102 is best explained by a 
pure H atmosphere mod el with T e g = 3830 K and logg 
= 8.08 (iKilic et al.ll20lJ l. These two ultracool WDs have 
masses of 0.77 and 0.62 Mq, respectively. Their total 
main-sequence + WD cooling ages and their kinematic 
properties indicate Galact ic halo membership . LHS 3250 
stands out in this sample. IBergeron fe Leggettl (l2002l l find 
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a best-fitting solution of T e ff = 3042 K, logg = 7.27, 
and log (H/He) = —2.7. Fig. [7] shows our model fits to 
the LHS 3250 SEP. Includi n g th e red wing of the Lya 
opacity dKowalski fe Saumonl 120061 1. we now derive slightly 
different parameters of T e g = 3064 K, logg = 7.38, and 
log (H/He) = —3.1. Even though the models predict deep 
absorption features around 0.8 and 1.1 /im, these features 
have never been observed in the actual spectra of cool and 
ultracool WDs. Clearly, the models have problems (likely 
due to problems with CIA calculations). Nevertheless, 
LHS 3250 is too bright compared to the pure H atmosphere 
models for ultracool WDs (see Fig. but its location in 
the colour-magnitude diagram and its SED are consistent 
with low surface gravity, mixed H/He atmosphere WDs. 

Fig. [5] shows our model fits to the J1102 SED us¬ 
ing pure H atmosphere models. We find T e ff = 3921 K 
and log g = 8.16 for this ultracool WD. These pa- 
ram eters are consistent with the analysis based on 
the Knwalskj_j£_S aumonl fl2006|) models within the errors 
dKihc_et_aL 201211. Both iKo walski fc Saumonl (l2006l i and 
iBergeron, Saumon fc Wesemaell ( 19951 ) models underpredict 
the flux in the i and 2 bands for this star. Addition of He 
provides a significantly better fit to the SED. The best¬ 
fitting mixed atmosphere model has trace amounts of helium 
with T e ff = 3327 K, log g = 7.64, and log (He/H) = —3.51. 
Based on this, all three of the previously analysed ultracool 
WDs with parallax measurements would be best explained 
by mixed H/He atmosphere models. 

The six ultracool WDs in our sample (see Fig. [ 6 ]) are all 
similar to LHS 3250. They are overluminous compared to 
pure H atmosphere WDs, and their SEDs and locations in 
colour-magnitude diagrams are matched fairly well by mixed 
H/He atmosphere WD models. The models have problems 
matching the peaks of the energy distributions and they 
predict absorption features at 0.8 and 1.1 /im that are not 
observed, but the overall fits are quite reasonable. Based on 
these, all nine (including the six in our sample) ultracool 
WDs with parallax measurements are best explained with 
H-rich (mixed) atmospheres. In addition, seven of the nine 
are low-mass objects with He-cores. These low-mass objects 
are about twice as large (see Table 4), and therefore four 
times as bright, as typical 0.6 Mg WDs. Hence, their over¬ 
abundance in the SDSS sample is not surprising. 


5.3 Common Proper Motion Pairs 

There are three common proper motion pairs in our sam¬ 
ple. These include two WD + WD pairs (J0747+2438 and 
J2239+0018) and one WD + K dwarf (J0045+1420). The 
lat ter was report ed as a common proper motion binary 
by iLuvtenl (ll987l) and iLepine fe Bongiorncl (120071 1. LSPM 
J0045+1421 (BD+13 99) is a G 8 V star 62.5 arcsec away 
from LSPM J0045+1420. We confirm that both the WD 
(J0045+1420) and the G 8 dwarf are at the same distance, 
making it a Sirius-like binary. Such binary systems can be 
used to constrain the initial-final mass relation. However, 
J0045+1420 is a 0.43 ± 0.07 Mq low mass WD with a cool¬ 
ing age of 2.7 Gyr. Due to its low-mass, J0045+1420 may 
itself be an unresolved binary, and it is impossible to con¬ 
strain its total age or the mass of its progenitor star. 

The remaining two WD + WD systems are also very 
useful as they provide a test of our cooling age estimates. 


J0747+2438 contains a 2.5lg;| Gyr old, 0.54 ± 0.05 M 0 DA 
WD and a 4.5/lQ'g Gyr old, 0.47 ± 0.05 M 0 DC WD that is 
best-fitted by a pure He atmosphere model. The separation 
between the two WDs is « 2000 au; it is safe to assume that 
both stars evolved independently. The cooling ages of the 
two stars differ by « 2cr, and this difference could be due 
to a difference in the mass of the progenitor main-sequence 
stars, though the lower mass star (J0747+2438N) is also 
the older star in this system. The progenitor stars of these 
relatively low-mass WDs were likely Sun-like stars that lived 
for 10 Gyr, and J0747+2438 is probably a very old binary 
system in the Galactic thick disc (it has a tangential velocity 
of only 38 ± 2 km s -1 ). 

J2239+0018 consists of a Teff = 4420 ± 90 K cool WD 
with a T e ff = 3510 ± 220 K ultracool WD companion. The 
two WDs are separated by only 1.85 arcsec, which corre¬ 
sponds to a physical separation of 155 au. Unfortunately, 
the relatively large error in our parallax measurement trans¬ 
lates into a large error in mass and cooling age estimates for 
this binary. J2239+0018B has a cooling age of G.S'Uj'l Gyr, 
whereas the ultracool WD J2239+0018A has a cooling age 
of 5.2lj '5 Gyr. Given the large errors, these two estimates 
are consistent within the errors. Further insight into under¬ 
standing similar binary systems will require more accurate 
distance measurements than those currently available. 


5.4 Disc versus halo 

Fig. [9] shows the mass versus temperature distribution of 
our parallax sample along with the theoretical isochrones for 
WDs with CO core composition s and thick envelopes, i .e., 
g(Hc) = 10 -2 and g(H) = 10~ 4 . [Bergeron, Leggett fe Ruiz! 
( 20011 ) explain the observed trend in the isochrones. The 
general trend is that low-mass WDs evolve (cool) faster 
than their counterparts, except that the onset of crys¬ 
tallization in the most massive WDs shortens the cool¬ 
ing times considerably. This leads to the parabola shaped 
isochrones. We also show the corresponding isochrones for 
the total main-sequence + WD cooling ages for r Js 
2 Gyr. We simply assume tMS = 1O(Mms //V / 0 )~ 2 ' 5 Gyr 
and Mms/M 0 = _81n[LMwp/M 0 )/O.4] (IWoodl Il992l ; 

iLeggett, Ruiz fc Bergeronlll998l) . 

The oldest globular cl usters in the halo a re currently 
producing 0.53 M 0 WDs dHansen et al , 2003). There are 
several WDs in our sample with M « 0.53 M 0 . J1436+4332 
is an excellent example of a potentially very old star. 
This DC WD is a T eff = 4630 ± 30 K, logg = 7.88, 
M = 0.50 ± 0.03 M 0 He-atmosphere WD with a WD cool¬ 
ing age of 5.4 ± 0.6 Gyr. Its progenitor main-sequence star 
was most likely a Sun-like star, with a main-sequence life¬ 
time of 10 Gyr. Hence, the total main-sequence + WD age 
of J1436+4332 could be as much as the age of the Universe. 

Fig. [9] demonstrates that WDs with M < 0.5 M 0 can¬ 
not have CO cores and also form through single star evolu¬ 
tion within the lifetime of the Universe. Hence, they must 
be either He-core WDs or unresolved double degenerate sys¬ 
tems. In fact, the majority of the targets in our parallax pro¬ 
gramme, including all of the ones with T e ff SC 4300 K, seem to 
be low-mass WDs with M < 0.5 M 0 . The isochrones shown 
in this figure are not appropriate for these low -mass WDs. 
Based on the lAlthaus, Serenelli fe Benvenuto! (l200ll ) mod¬ 
els, they have cooling ages ranging from 2.9 to 9.7 Gyr. Since 
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Figure 9. Masses of all stars in our sample as a function of effective temperature. The black dots correspond to the WDs with fi¬ 
nch atmospheres while the white dots represent the He-rich WDs. The red dots represent the 15 cool and ultracool WDs with mixed 
atmospheres. Also shown are theoretical isochrones with their ages labelled in Gyr; solid lines correspond to WD cooling ages only, while 
the dotted lines also include the main-sequence lifetime. 


the prior evolution of these systems, including the masses 
of their main-sequence progenitors, is unknown, their total 
ages cannot be estimated. However, some of the low-mass 
ultracool WDs are clearly very old. 

Tables [2] and [4] and Fig. [TO] present tangential veloci¬ 
ties and cooling ages for our parallax sample. The cooling 
ages range from about 2 to 10 Gyr, whereas the tangen¬ 
tial velocities of all but two of the targets are less than 
150 km s -1 . The only targets that display halo kinemat¬ 
ics are J0301—0044 and J1107+4855 (hereafter J0301 and 
J1107, respectively). J0301 and J1107 have tangential ve¬ 
locities of 167-192 km s“\ T eff « 4500 K, M = 0.45- 
0.50 M©and WD cooling ages of 5-6 Gyr. If both stars are 
single CO core WDs, their progenitors would be Sun-like 
stars with main-sequence lifetimes of 10 Gyr. Hence, their 
total main-sequence + WD cooling ages would be ~14 Gyr, 
which is consistent with a halo origin. 

In Table[2j the last three columns list the (U, V, W) com¬ 
ponents of the space velocities for each WD. These space 
velocities have been computed by combining the observed 
parallaxes and proper motions for each W D using the pre¬ 
scription of I Johnson fc Soderbloml d 19871 '). Since we do not 
have any radial velocity measurements for our WDs, we as¬ 
sume a radial velocity u ra d = 0 km s -1 in the calculation of 
U. In Fig. 1111 we plot the resulting W versus V (top), and 
U versus V (bottom) velocity distributions. We also include 
the 2cr contours for the Galactic thin disc (do tted), thick disc 
(dashed) and stellar halo populations (solid) dChiba fc Beers! 


l2000l ). It is clear from Figure fill that both J0301 and J1107 
are most likely halo WDs. The distribution of the remain¬ 
ing sample, including the IR-faint (the DC WDs with mixed 
H/He atmospheres) and ultracool WDs, is consistent with 
disc membership. The ratio of thick to thin disc WDs is 
18/34. If we assume that J0301 and J1107 are indeed bona 
fide members of the halo, the observed velocity distribution 
suggests 63 per cent/33 per cent/4 per cent proportions for 
the contribution of thin disc, thick disc and halo WDs for 
our sample. 

We note that this is the first time a large number of 
ultracool WDs have distance and tangential velocities avail¬ 
able, and in contrast to the expectations, they all seem to 
be members of the disc. These WDs provide an independent 
constraint on the thick disc population; the Galactic thick 
disc is at least 10 Gyr old. 


6 CONCLUSIONS 

We present parallax observations of 54 cool and ultracool 
WDs. Our sample includes one new WD within the local 
25 pc sample and five stars within 30 pc. All but two of 
them have tangential velocities smaller than 150 km s^ 1 . 
J0301—0044 and J1107+4855 are the only objects in our 
sample with kinematics and ages that are consistent with 
halo WDs. The rest of the objects, including the ultracool 
WDs, are members of the Galactic thick disc. The oldest 
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Figure 10. The tangential velocity, utan, plotted as a function 
of the WD cooling age. The black dots correspond to the WDs 
with H-rich atmospheres while the white dots represent the He- 
rich WDs. The red dots represent the 15 cool and ultracool WDs 
with mixed atmospheres. 
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WDs in this sample have WD cooling ages of 10 Gyr, pro¬ 
viding a firm lower limit to the age of the thick disc. Many 
of our targets are low-mass WDs. These are either single 
He-core WDs or unresolved double degenerates. It appears 
that we have detected the brighter population of cool and 
ultracool WDs in the solar neighbourhood, and the fainter, 
normal CO core ultracool WDs remain to be discovered in 
large numbers. Future and upcoming astrometric surveys 
such as the LSST will find those fainter and more massive 
ultracool WDs. 
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Figure 11. W versus V (top) and U versus V (bottom) veloc¬ 
ity distributions for the 54 WDs in our sample. The ellipsoids 
denote the 2<r contours for the Galactic thin-disc (dotted), thick- 
disc (dashed) and stellar halo populations (solid). The black dots 
correspond to the WDs with H-rich atmospheres while the white 
dots represent the He-rich WDs. The red dots represent the 15 
cool and ultracool WDs with mixed atmospheres. 
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